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Abstract 

We report the results of a Monte Carlo study of the continuum 
limit of the two dimensional 0(3) non-linear a model. The notable 
finding is that it agrees very well with both the prediction inspired by 
Zamolodchikovs’ S'-matrix ansatz and with the continuum limit of the 
dodecahedron spin model. The latter finding renders the existence of 
asymptotic freedom in the 0(3) model rather unlikely. 

Asymptotic freedom is one of the corner stones of modern particle physics. 
This special property of certain non-Abelian models is supposed to explain 
many remarkable properties of quantum field theory, such as the existence 
of a non-trivial continuum limit, the possibility of grand unification of all 
interactions, etc. This property was discovered in 1973 in perturbation the¬ 
ory (PT) [|, 1. However the correctness of this technique within a non- 
perturbative setting, such as the one offered by lattice quantuum field theory, 
remains mathematically unsettled. In the past, we have repeatedly pointed 
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out that there were good reasons to doubt the correctness of the PT approach 
in non-Abelian models [!3|, [|. |j]. 

In this letter we would like to report some new numerical results which, in 
our opinion, make the existence of asymptotic freedom in the two dimensional 
(2D) non-linear 0(3) a model extremely unlikely. Namely we find that the 
massive continuum limit of this model is, as far as we can see, identical to 
that of the dodecahedron spin model, which as a discrete spin model must 
undergo a transition to a phase with long range order (Iro) at a finite inverse 
temperature (3 and hence it cannot possibly be asymptotically free. The 
only escape out of this bind would be the following scenario: the presently 
observed agreement between the 0(3) and the dodecahedron spin models is 
accidental and the phase transition in the dodecahedron model is actually 
weakly first order, a possibility which we cannot rule out, yet for which there 
is no theoretical basis. 

The 0(3) model is believed to possess only one phase, with exponentially 
decaying correlations. It is of interest because besides being believed to 
be asymptotically free, it is a model for which several theoretical predictions 
have been made. They stem from a conjecture put forward by Zamolodchikov 
and Zamolodchikov [6] regarding its S'-matrix. Two predictions were made 
using this S'-matrix ansatz: 

- Hasenfratz, Maggiore and Niedermayer (HMN) [7] used it in conjunction 
with the thermodynamic Bethe ansatz to predict the ratio A/m 

- Balog, using the ideas of Karowski and Weisz about constructing form 
factors out of the Zamolodchikov S'-matrix || , extended the work of Kirillov 
and Smirnov || and computed the current 2-point function and then 


later, together with Niedermaier and Hauer H. also the 2-point function of 
the energy-momentum tensor as functions of p/m. Balog and Niedermaier 
jT2i extended this construction to the two-point functions of the spin and the 
topological density. 

While it has been known for several years that the HMN prediction for 
A/m disagreed with the Monte Carlo value by at least 15%, as far as we 
know no previous work exists comparing the predictions for the current and 
spin 2-point functions with Monte Carlo data. In this letter we report the 
results of some Monte Carlo investigations of these issues. For p/m less than 
about 20 we find excellent agreement with the theoretical predictions for 
the current and spin 2-point functions. On the other hand, even though we 
reach correlation length £ = 167 the value of A/m disagrees with the HMN 
prediction by 15% and displays an increase with £, suggesting that it will not 


2 





converge to the predicted value. 

The numerics were performed using a variant of Wolff’s one cluster al¬ 


gorithm (T^j. One run consisted of at least 100,000 cluster updates used for 
thermalization, followed by at least 7 bins of 100,000 clusters used for taking 
measurements. Thus each run consisted of up to 1,000,000 clusters used for 
measurements. This procedure was repeated, starting with a different initial 
configuration; the total number of measurement clusters is recorded in Tab.l. 
The error was estimated in each run from the bins. If more than one run 
was performed, we also estimated an error by treating each run as a bin. 
Typically the latter error was larger and was recorded as the actual error 
estimate. 

The quantities measured were: 

- Spin 2-point function G(p ): 


g ( p ) = ^( \m 2 y, 

Current 2-point function J(p): 


s(p) = y^e ipx s(x) 


( 1 ) 


J(p) = T2 £ £<l i»r>; 3» = £e^>) 


( 2 ) 


a= 1 (i=l 


where j*{x) = (3e abc s b (x)s c (x + fi) 
- Magnetic susceptibility y: 


X = G(0) 


Correlation length £: 

e = 


2 sin(7r/L) 


\J(G(0)/G(1) - 1) 


( 3 ) 


( 4 ) 


We measured these observables in both the 0(3) and dodecahedron spin 
model. However only in the 0(3) model there exists a continuous symmetry, 
leading to a conservation law and to Ward identities Jl4| . The normalization 
of the current in eq.(2) was fixed by the latter. In the dodecahedron spin 
model there is no a priori conserved current nor any Ward identities and it 
is a remarkable finding of our investigation that in the continuum limit the 
current 2-point function in the dodecahedron spin model is identical, up to 
an overall normalization, with that in the 0(3) model. 
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Let us briefly review what is known about the dodecahedron spin model 
with standard nearest neighbour interaction: 

a) It possesses a high temperature phase with exponential decay of correla¬ 
tion functions. 

b) It possesses a low temperature phase with Iro and at least one pure phase 
for each dodecahedron vertex. 

c) It seems to have an intermediate phase with algebraic decay of correlation 
functions. 

While the first two properties follow easily from convergent high and low 
temperature expansions, the last one has not been proved rigourously. It 
would follow from a rigorous inequality derived by Richard and us [ 151, con¬ 


necting the dodecahedron model with the Z{ 10) model, provided one knew 
that 

Pm(D) > (3 C {Z{ 10)). (5) 

Here (3 m (D) is the inverse temperature below which the dodecahedron model 
exhibits Iro and (3 C (Z( 10)) the inverse temperature above which Z(10) ex¬ 
hibits exponential decay. In ref. |I6| we gave numerical evidence for the 


existence of such an intermediate massless phase. In such a phase, one could 
expect the discrete icosahedral symmetry to be enhanced to 0(3), just as 
Z(N) for N > 4 is enhanced to 0(2) [|T7|. The results reported in this pa¬ 
per regarding the similarity of the massive continuum limits of the 0(3) and 
dodecahedron model suggest that indeed this is the case. 

We investigated the massive continuum limit in a thermodynamic volume. 
To be thermodynamic one needs L/£ >7 and all of our runs obyed this 
criterion. To reach the massive continuum limit one takes £ —> oo keeping 
p£ fixed. Since numerically one cannot increase £ to oo, one must control 
lattice artefacts. We did this in two ways: 

1. For a given action, we studied the approach to the continuum limit at 
fixed p/m as a function of £ (m = l/£). 

2. We studied the dependence of the results upon the lattice action. 

For the standard nearest neighbour action we varied the correlation length 
from 11 to 167. Besides this action, we took measurements also with the 


‘cut-action’ (/? = 0 but grad(s) restricted |18|) and with a ‘Gaussian perfect 
action’ ( s{i) ■ s(j)) gets replaced by the the square of the angle between the 
two spins and there is coupling to the nearest neighbour and to the site 
diagonally located |L9|). Our main conclusion is that the continuum limit is 
approached more slowly at increased p/m and that for instance at p/m =50 
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one needs £ around 100 for a 1% deviation in J(p). We observed no striking 
effect related to the lattice action used, only perhaps a slight advantage 
gained with the ‘Gaussian perfect action’ for the dodecahedron spin model 
(for which the Gaussian improvement is totally ad hoc). 

To motivate the presentation of our results let us briefly discuss the form- 
factor approach of Balog et al. In the spin or current 2-point function they 
insert the complete set of n-particle states and write them as an infinite sum, 
given in terms of form-factors. O(N) invariance dictates that in G(p) only 
odd terms contribute, while in J(p) only even ones do. Thus G(p) starts with 
the 1-particlc contribution, which does not depend on the S'-matrix. Thus, 
to see the dynamics, it is better to subtract this contribution; so we report 

G s {x) = [G(x)/G{0) - l/(x 2 + 1)/1.001687] x x 2 (6) 


Here x = p/m , the division by G(0) represents the wave function renormal¬ 
ization necessary to have a nontrivial continuum limit and the numerical 
factor 1.001687 represents the difference between our and Balog’s normal¬ 
ization of the spin 2-point function [20]. The whole expression was mul¬ 


tiplied by (p/m) 2 because, according to renormalized perturbation theory, 
after such multiplication, the expression should diverge as log(p/m). The 
current 2-point function requires no renormalization (in 0(3)) and, if one 
accepts renormalized perturbation theory, is expected to diverge as log(p/m) 
at large p |TD[ . In the dodecahedron spin model, the normalization of the cur¬ 
rent is arbitrary, hence we could either consider J(p/m)/J( 1) and compare it 
with the same quantity measured in 0(3) or simply see if after multiplication 
by a constant J(p) in the two models coincide. 

Results: 

Perturbative renormalization group arguments lead to the predictions 


£ 


C £ e 


27T/3 


2tt/3 


-(! + £ 


k> 1 




x 


C x e 4 *% 

(27T/5) 4 


£ 

k> 1 


/ 1 * 


(7) 


the action dependent coefficients a*, and bk are known for the standard action 
up to 4 loops [22|. The HMN prediction is that 


C £ = -2~ 5/2 e-” /2 « 0.012487 
? 8 


( 8 ) 


Acoording to ref. |2^. |2T[ the HMN prediction implies also C x = 0.0145. 
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Figure 1: C^ = A/m compared to the HMN prediction (horizontal line). The 
data points correspond upwards to 2 loop, 3 loop and 4 loop PT, respectively. 


In Fig.|T] we present our determination of A/m in the 0(3) model with 
standard action, in Fig.|| our determination of O x , compared to the predic¬ 
tions by [i7| and 0 , respectively. The values measured and the number of 
clusters used are recorded in Tab.l. 

In Fig.[i] and |3]) we compare G s (p/m ) (eq.(6)) and J(p/m) (eq.(2)), re¬ 
spectively, for the 0(3) model with standard action at different correlation 
lengths with the prediction of Balog and Niedermaier |T2| and Balog |T0| . 
respectively. The errors are less than 1CT 3 and therefore not shown. 

Tab.l: Monte Carlo data for the standard action 


p 

L 

clusters 

£ 

X 

1.5 

78 

1 x 10 6 

11.03(2) 

175.75(43) 

1.6 

140 

1 x 10 y 

18.97(4) 

447.30(6) 

1.7 

250 

2 x 10 6 

34.55(7) 

1268.(3.) 

1.8 

500 

4 x 10 6 

64.67(8) 

3831.(6.) 

1.9 

910 

2 x 10 6 

122.15(39) 

11847.(67.) 

1.95 

1200 

7 x 10 5 

168.39(71) 

21148.(123.) 


In Fig.|| (|6) we compare G s {jp/m) ( J(p/m )) for the 0(3) and the dodec- 
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Figure 2: C x compared to the prediction of Alles et al (horizontal line). The 
data points correspond upwards to 2 loop and 3 loop (indistinguishable) and 
4 loop PT, respectively. 


ahedron models, both with standard action. For the dodecahedron model 
J(p/m ) was multiplied by a renormalization factor determined by fitting. 

Discussion: 

1. The data in Fig.l agree very well with prior data |24|, ^5f and indicate 
that most likely A/m will increase past the HMN prediction. Attempts to 
improve this situation by using the so-called ‘energy scheme 1 |23| produced 
similar results (the curve is shifted closer to the predicted value but it shows 
a non-vanishing slope away from it). 

2. A comparison of \(f3) with the predictions of asymptotic freedom shows 
the same behavior: the data increase faster than expected, in agreement with 
the existence of a critical point at finite /3. From the data in Tab.l we find 
that C x is growing steadily, so far reaching 0.0111, compared to the ‘exact’ 
prediction of a constant value of 0.0145. 

3. The numerical agreemeent of G s {p/m) and (the normalization ad¬ 
justed) J(p/m ) in the 0(3) and dodecahedron spin model and their excellent 
agreement with the Balog and Balog-Niedermaier prediction strongly suggest 
that the two models possess the same continuum limit, which is correctly pre- 


7 







B 

\ 

o" 


p/m 

Figure 3: G s (p) for 0(3) as a function of p/m. The solid line is the prediction 
of Balog and Niedermaier. 



dieted by the Zamoloclchikovs A-matrix ansatz. Since the massive phase of 
the dodecahedron model terminates at finite (3, its continuum limit cannot 
have asymptotic freedom. Indeed, for instance the running coupling g r (L) 
introduced in |26| cannot vanish at short distances, because at any [3: (a) 


g r {L) 0 (2) does not vanish and (b) g r (L)o( 3 ) > g r (L)o{ 2 ) 

4. Balog and Niedermaier ||28|| have put forth an ansatz for the contribu¬ 
tion of higher n-particles to G{p/m) and J(p/m). According to them this 
ansatz would explain why the contribution of n > 6 is negligible at p/m < 10 4 
and how the same contributions sum up at asymptotic p/m to produce the 
logarithmic increase of both G s (p/m) and J(p/m ) corresponding to asymp¬ 
totic freedom. Our data suggest that their ansatz is inadequate, since the 
agreement of 0(3) with the dodecahedron model rules out such a logarithmic 
increase in J{p) 


see 
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Figure 4: J(p) for 0(3) as a function of p/m. The solid line is the prediction 
of Balog and Niedermaier. 
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